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The stickiness effect suffered by chaotic orbits diffusing in the phase space of a dynamical system is studied in this paper. Previous 
works have shown that the hyperbolic structures in the phase space play an essential role in causing the stickiness effect. We 
present in this paper the relationship between the stickiness effect and the geometric property of hyperbolic structures. Using a 
two-dimensional area-preserving twist mapping as the model, we develop the numerical algorithms for computing the positions of 
the hyperbolic periodic orbits and for calculating the angle between the stable and unstable manifolds of the hyperbolic periodic 
orbit. We show how the stickiness effect and the orbital diffusion speed are related to the angle. 
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1 Introduction 

In a 2n-dimensional phase space of a dynamical system, 
given long enough time, a chaotic orbit would visit any 
place in the connected chaotic region, owning to the er- 
godicity property. Only those regions surrounded by the 
(2n - l)-dimensional invariant surface (curve) are forbidden 
to the chaotic orbit. During its diffusion in the phase space, 
a chaotic orbit generally spends much longer time wandering 
near the regular regions than elsewhere. This phenomenon is 
called the "stickiness effect" [1,2], and it causes the anoma- 
lous chaotic transport [3-7]. The phase spaces of most phys- 
ically realistic systems consist of chaotic and regular com- 
ponents, thus the stickiness effect is a generic phenomenon. 
Especially, in those nearly integrable systems such as the So- 
lar System [8], the slow diffusion is a key point to understand 
the stabilities of those systems. As a matter of fact, the dif- 
fusion (transport) in the phase space is a very important issue 



in many areas, with interests not only of the general theory of 
dynamical system but also of direct applications in physics 
and astrophysics. For example, this type of dynamics has im- 
plications in the tokamak fusion [e.g. 9], turbulence [e.g. 10], 
plasma physics [e.g. 11], semiconductor superlattices [e.g. 
12], quantum chaos [e.g. 13], galaxy dynamics [e.g. 14], etc. 

The stickiness effect was first recognized in the close vicin- 
ity of the KAM tori [1, 2], but it seemingly may arise from 
many other structures in the phase space, including the is- 
land chains, the cantori, the hyperbolic structures, and the 
asymptotic curves of the unstable periodic orbits [see e.g. 2, 
15-24]. We have shown in our previous papers [25, 26] that 
the hyperbolic invariant sets play a critical role in causing the 
stickiness effects. Using models of both 2-dimensional (2D) 
mapping [25] and 3-dimensional (3D) mapping [26], we have 
shown that a diffusing orbit is generally "stuck" by the hyper- 
bolic structures in the 2D or 3D phase space. There seems to 
be an essential relation between the hyperbolic structures and 
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the stickiness effects. 

Many studies on the stickiness effect focused on the diffu- 
sion behaviors in the phase space, by illustrating what struc- 
tures may have stickiness effects [see e.g. 16, 20, 21, 23], or 
by describing the rules of diffusion in different regimes [e.g. 
8, 16, 17, 27]. As for the geometric details of the structures 
causing the stickiness effect, Contopoulos and Harsoula [23] 
calculated the sizes of gaps in a cantorus being crossed by 
an orbit and relate the size of the "maximum gap" with the 
diffusion time. 

Obviously the diffusion orbits have to pass through the 
"gaps" when crossing a cantorus. In a 2D phase space, the 
hyperbolic periodic orbits are included in cantori and island 
chains. When an orbit crosses these structures on its diffusion 
route, it approaches from one side to the hyperbolic periodic 
point along the stable manifold and then it may leave from 
the other side along the unstable manifold. Therefore, the 
geometric features of the phase space influence strongly the 
characteristic of diffusing trajectories. In this paper, we try 
to explore the geometric structures around the cantori, in par- 
ticular, we study the angle between the stable and unstable 
manifolds around the cantori. 

To do this, it is necessary to know first the exact positions 
of the cantori and the directions of the stable and unstable 
manifolds mentioned above. But, since a cantorus is unsta- 
ble, it's difficult to locate its position precisely. Consequently, 
the directions of the manifolds can hardly be determined ac- 
curately. Using a 2D mapping model, we will introduce the 
methods for calculating such positions and directions with 
high precision, then calculate the angle between the direc- 
tions of the stable and unstable manifolds, and finally we will 
study the relation between the diffusion speed and the sizes 
of the angles. We find that the diffusion speed increases when 
the size of angles increases. 

The mapping model used in this paper is obtained by a 
modification to the so called "standard mapping". Contrary 
to the standard mapping whose "twist" depends linearly on 
the action-variable, our mapping model has a nonlinear twist. 
Such a nonlinearity in twist is often seen in physics [28]. 

The paper is organized as follows. In Section 2, we intro- 
duce the mapping model. We show in Section 3 how to lo- 
cate the precise positions of hyperbolic periodic orbits in the 
phase space and present the method of calculating the angle 
between stable and unstable manifolds. The characteristic an- 
gles of the hyperbolic structures in different regions and the 
diffusion speeds of orbits diffusing across these regions are 
compared in Section 4. And finally we make the conclusions 
in Section 5. 

2 Mapping Model 

Generally a mapping is much easier to handle than a set of 
differential equations (e.g. a Hamiltonian system) as a dy- 
namical system. Thus we use a mapping model in this paper 
as usual. 
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Figure 1 The phase space of the mapping (Eq. (6)) with the parameter 
k = 20.75. Some hyperbolic periodic orbits (fixed points) are located and 
plotted. To present a clean appearance, the period-two orbits are indicated 
by crosses, and the period-3, 4 orbits are indicated by triangles and squares, 
respectively. The 3/1 and the 35/11 periodic orbits are indicated by a solid 
dot and stars. The rotation number of the orbits are labeled on the right side. 



2.1 The mapping 

We begin from the well-known standard mapping [29]: 



x — x — y , 
5' = y + h sin(27r£), 

which can be generated from a generating function 
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(x-x) 
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cos(27ti), 



through equations 



y = Sz = m = x-x'-£; sin(27tx), 
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(3) 



Here we modify the generating function by adding an extra 
high-order term Ux — x') 4 , 



Ilk 
S = -(x - x'f + -(x - x') A + — - cos(2ttjc), 
2 4 An 1 



(4) 



thus the corresponding area-preserving mapping is given by 

(5) 

mod(l), 



y = S x = (x - x 1 ) + (x - x') 3 - sin(27t^), 
y = S x , = (x-x') + (x-x') 3 . n 

Explicitly, the mapping reads 

( x> = x-^ + ^t-^-^l 

I V =y+h sin(27Dc). 

(6) 

This mapping is defined on the cylinder < x < 1, -co < 
y < +oo, and k is the only one perturbation parameter. As 
in the standard mapping, the y in this mapping is the action 
variable, and the x is the angle variable that can be performed 
the modulus over 1. As usual, in this paper, by "diffusion" 
we mean the drifting of the action variable y in the phase 
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space. In the study by Cheng and Sun [30], it was shown 
that for any given k, this mapping behaves chaotically around 
\y\ < M (where M > 0), and the KAM curves exist at large |y|. 
When k is small, the phase space of the mapping is mainly 
full of (horizontal) invariant tori. As k increases, the KAM 
tori break, first at low \y\ value, and chaos sets in. The area 
occupied mainly by chaotic orbits extends outward from the 
region around y = 0. We show in Fig.Q~]the phase space when 
k - 20.75. Note that the negative y half of the phase space is 
symmetric to the positive half that has been shown in Fig.Q] 
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Figure 2 The rotation number of the mapping (k = 20.75). Open squares 
show the change of rotation number along the vertical line x = 0.0, while the 
crosses are for the line x = 0.5. 



The mapping is an area-preserving monotone twist map- 
ping, thus along a vertical line in the phase space the rotation 
number increases monotonically. The rotation numbers along 
two vertical lines are computed and shown in Fig. [2] Clearly, 
the scattering distribution in the lower left corner of Fig. [2]is 
a reflection of the chaotic motion at low y value in the phase 
space. The "plateaus" in both cases arise from the fact that 
the vertical lines cross stable islands, on which the rotation 
numbers are constants. 



2.2 The orbital diffusion in the phase space 

Any chaotic orbit in the connected chaotic sea of the phase 
space will wander ergodically as time tends to infinity. But, it 
is known that an orbit spends much longer time in some spe- 
cific regions than elsewhere, e.g., the close vicinity around an 
embedded island, the region where a KAM torus had newly 
broken, the region occupied by hyperbolic structures, etc. 
This is called the stickiness effect, as we mentioned above. 

To show a typical orbital diffusion process and the stick- 
iness effect, we arbitrarily select an initial point (xo,yo) = 
(0.49, 29.0) not far away from the regular region in the upper 
part of Fig. Q] and follow its evolution. The diffusion of this 
orbit is shown in Fig. [3] 
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Figure 3 The diffusion process of an arbitrary orbit starting from (xo,yo) = 
(0.49, 29.0). The red, black and blue points represent its trajectories during 
the time (iteration number n) intervals of (1 , 9 X 10 4 ), (9 X 10 4 , 4.5 X 10 7 ) and 
(4.5 X 10 7 , 7 X 10 7 ). After ~ 6 X 10 9 iterations, it finally crosses the compact 
partial barrier (see text), which makes the bottom boundary of blue points, 
and enters the open chaotic sea around y = (green dots). 

On its route of downward diffusion, the orbit meets some 
structures that play the role of obstacles. As shown in Fig. [3] 
the orbit wanders around the 3/1 fixed point and the corre- 
sponding island (see Fig. Q] to find the positions of island- 
chains of different rotation numbers) at beginning (red points 
trajectory). This island does not show a significant stickiness 
effect, but the orbit spends much longer time in next stage 
around the vicinity of the 11/4 island-chain (black dots in 
Fig. [3). After that, it was obstructed by another barrier, which 
prevents the orbit from further diffusing downward, and only 
after pretty long time, the orbit can reach the 8/3 island chain 
(blue dots). Then, the diffusion was strongly obstructed by a 
partial barrier below the 8/3 island chain. Only after ~ 6x 10 9 
iterations, the orbit can cross this "compact" barrier and en- 
ter finally the open chaotic sea around y = (green dots). In 
fact, the rotation number around this compact partial barrier 
is around 2.6 (see Fig. [2}. We suspect that this difficult-to- 
cross barrier is the remnant of the last and most robust KAM 
torus that has the "golden rate" rotation number. And this 
will be verified in next section. 

One thing we would like to address here is that the or- 
bit may cross an obstacle in both forward and backward di- 
rections, that is, an orbit may cross the same partial barrier 
several times back and forth. After leaving a region, it may 
re-enter the region again some time later. That's why we see 
the black dots and red dots are mixed around the 3/1 island in 
Fig. [3] 

Typically, an orbit starting above the lowest and difficult- 
to-cross barrier mentioned above would spend time as long as 
several 10 9 iterations before crossing it. In attempt to show 
some details of this barrier, we plot in Fig. |4]three orbits start- 
ing close to it, from initial points Pi = (0.509, 21.667), P2 = 
(0.509, 21.672) and P 3 = (0.509, 21.677) respectively. These 
initial points are indicated by three crosses in Fig. |4j). From 
this magnified picture, we realize that P3 (black) is on an 
island-chain with a high period number. The orbit initialized 
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at P2 (red) finally diffuses upward after several 10 9 iterations, 
while the one from Pi (blue) diffuses downward. There is a 
mixed region where we can find both red and blue dots, in- 
dicating clearly the existence of the difficult-to-cross partial 
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Figure 4 The most sticky region. Three initial points (crosses) and the 
corresponding trajectories (dots) are shown in black, red and blue. 



3 Angle between stable and unstable manifolds 

The diffusion of an orbit is through the hyperbolic structures 
between the islands in an island-chain, which is in fact the 
relic of the KAM tori if the rotation number in this position 
does not satisfy the Diophantine condition [31]. On the other 
hand, if the rotation number is irrational enough, there would 
be a cantorus left after the breaking of the KAM torus. And 
an orbit's diffusion is through the gaps in the cantorus. But 
the cantorus itself is composed of hyperbolic points and in 
the gaps exist the hyperbolic structures. 

In this paper, we will describe the geometric property of 
the hyperbolic structures and find out how this property af- 
fects the diffusion speed. We will show that the diffusion 
speed is determined by the angle between the stable and un- 
stable manifolds of the hyperbolic structure. In this section, 
we introduce the algorithm for calculating the angle between 
the stable and unstable manifolds of a hyperbolic structure. 



3.1 The locations of hyperbolic periodic points 

Usually a recursive algorithm is employed to compute the lo- 
cations of periodic orbits. But it does not work for the hyper- 
bolic ones because these orbits are unstable. Fortunately, the 
hyperbolic periodic orbits we are looking for in this model 
are of minimal action, therefore the gradient algorithm could 
be used to find these orbits. The existence of these minimal 
periodic orbits is guaranteed by Aubry-Mather theory [32]. 

Suppose xo, x\, ■ ■ ■ ,Xp-\,x p are the variables in the con- 
figuration space of a minimal periodic orbit of rotation num- 
ber q/p, with q and p being coprime integers. Obviously, 
x p = xo + q. The action functional / corresponding to this 
periodic orbit is: 

/(x , xi, ■ ■ ■ , Xp-\) = S (x , xi) + S (xi, xz) + ■ ■ ■ + S (x p -\ , x p ), 

(7) 

where S is the generating function (Eq. (4)). Then the min- 
imal periodic configuration (xo, xi, ■ • ■ , x p -i,x p ) is the mini- 
mum of /. Therefore, a periodic orbit of rotation number q/ p 
can be found by locating the minimum of the above men- 
tioned functional. To find the minimum, we may use the gra- 
dient algorithm. 

Consider the following differential equation: 



a = -v/(H) 



(8) 



where H = (ft, ft, ■ ■ ■ iftj-i) r is a column vector and the dot 
over it indicates time derivative. A trajectory defined by this 
differential equation converges to a minimum of the func- 
tional / as t — > oo. Thus, a minimal periodic orbit can be 
obtained through integrating this differential equation. 
Specifically in our model, since 

S Xk (x k -i,x k ) = -y k , S Xt (x k ,x k+ i) =y k ,k=\,2,...,p. (9) 

where the subscript denotes the partial derivative over the cor- 
responding variable, and 

fx t = S xi (x k -\,x k ) +S Xt (x k , Xfc+i) = -y k + y k = 0, (10) 

the corresponding differential equations reads: 



ft 



-(ft- 

-(ft- 



ft) - (ft - ft) 3 + h sin(27^ ) + [ft-i - (ft + ?)] + [ft- 
ft) - (ft - ft) 3 + |f sin(27t^) + + (ft - £) 3 , 



(ft + <7)] 3 



ft = -(ft - ft + i) - (ft - ft +1 ) 3 + j- sin(27Tft) + (ft_! - ft) + (ft_i - ft) 3 



(11) 



ft- 



"(ft-2 
"[ft-i 



■ft-l) - (ft-2 - 
■ (ft + <?)] - [ft- 



ft-i) 3 + -t sin(27Tft_ 2 ) + (ft_ 3 - ft 
i - (ft + <z)] 3 + h sin(27tft_i) + (ft 



-2) + (ft-3 
-2- ft-l) H 



"ft 

(ft 



-2)\ 
-2 - ft- 



) 3 . 



The solution of these differential equations can be computed 
numerically. We adopt the 7th order Runge-Kutta-Fehlberg 
algorithm (RKF7(8)) to integrate the differential equations to 
time 10 4 , unless the solution reaches the convergent value 
within a given error when we regard the solution attains the 
required accuracy. In this paper, we set the controlling error 
to be 10~ 14 , that is, if the difference between the successive 



values obtained from the integrator is smaller than 10 14 , we 
regard the convergent solutions have been achieved. 
In fact, the solution converges quite quickly. We show in 
Fig.|5]an example of the convergence of the trajectory to the 
minimal periodic orbit of rotation number 9/4. We start from 
(ft = 0.05, ft = 0.30, ft = 0.55, ft = 0.80), and the so- 
lution converges to the final value very soon, at integrating 
time f ~ 11. 
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Figure 5 The convergence of the solution to the periodic orbit of rotation 
number 9/4. 

With the x values of the minimal periodic orbit, the cor- 
responding y values can be computed using the generating 
function, via Eq. (5). And finally we get the positions of the 
minimal periodic orbit, i.e. the hyperbolic periodic orbits in 
the phase space. We have over-plotted some of these calcu- 
lated points on the phase space as shown in Fig.Q] 



3.2 Characteristic angle of the hyperbolic structure 

Knowing the positions of hyperbolic periodic orbits, we now 
turn to their properties. One of the most important characters 
of the hyperbolic structure is the angle making by directions 
of the stable and unstable manifolds, which will be called 
"characteristic angle" for short hereafter. The directions of 
stable and unstable manifolds are determined by two special 
tangent vector fields along the periodic orbit, while the tan- 
gent vector fields can be obtained by a limit processes. Along 
an orbit, the projection of the tangent vector field onto the 
horizontal direction gives a Jacobi field, which can be easily 
calculated by using the generating function of the mapping. 
We start our computations with the definition of Jocobi fields: 
Assume a hyperbolic periodic orbit {x^ = _ x of rotation num- 
ber q/p has been found, and it satisfies 

xo, x\, ■ ■ ■ ,x p -\ and x n+p — x„ + q (n = 0, ±1, ±2, • • •). 

Then, {ft}," is said to be a "Jacobi field" along the orbit 



S 12(*n-l,X„)ft_i + [S22(x n -\,X„) + S u(X n , X„+{j\ ft + S I2(x„, X„+l)ft + i - 0, 



(12) 



where S is the generating function Eq. (4) and the subscript '1' and '2' indicate the partial derivatives with respect to the first 
and second variable of the function. Actually, the above equations can be written in a matrix form. Set 



A„ = S 22(x„-i, x n ) + S u (xn, x„+i), 
Then we may write the equations as: 



B n — S nix,,, x n+ \), 77 = 0, +l,+2, ■ 



(13) 



B„_i An B„ 



ft-i 
ft 
ft+i 



(14) 



Any finite piece of this triangular matrix is positive definite 
if the corresponding configuration \x n ] is minimal. A Jacobi 
field along a minimal configuration induces an orbit of the 
corresponding tangent map [33]. 

To compute the Jacobi fields along a minimal periodic orbit 
{x n }, which correspond to the stable or unstable directions, we 
can get the approximation of the real Jacobi fields by solving 
the following linear equations: 



( A] B\ 
Bi A 2 

6 ••• 



o 

B 2 



B„-2 





A„-\ 
B„-\ 







B„-\ 
A„ 



( 6 1 




( - fl o 'I 


ft 







ft-1 







v ft > 




I ) 



(15) 

As « — » oo, we get the real Jacobi field. Note that £i in above 
equations is monotonically increasing as « goes to infinity. In 
fact, the limit value of %\ is the projection on x-axis of the 
stable direction at point {x\,y\). Since (from Eq. (5)) 

yo = Si(x ,x{), 

the projection on y-axis of the stable direction at point 
(xo, yd) can be computed from 



with £o = 1- Finally the stable direction at (xo,yo) is de- 
fined by i s = (1,770). Note that this Jacobi field could be 
used to calculate the Lyapunov exponent A of the orbit, where 




Figure 6 The convergence of fi in Eq. (15). This is the case for the 
minimal periodic orbit of rotation number 9/4. (a) ^\ versus recursive num- 
ber n; (b) The error versus n, where the error is defined as the difference 
Afi=6(n)-fi(1000). 

We calculate ft (77) using a recursive algorithm. An exam- 
ple for the minimal periodic orbit of rotation number 9/4 is 
shown in Fig. [6] The convergence of £1 with the increasing 
recursive number n can be seen clearly. 
In a similar way, reversing the direction of calculation, we 
obtain which is the projection of unstable direction on 
x-axis at point (x_i,y_i). Considering 

yo = -S 2 (x-i,xo), 

we compute the y-axis projection of the unstable direction at 
point (xo,yo) through: 



770 = S n(x , xi)£ + S u(x , xi)£i, (16) 



n'o = Sv2(x-i,Xo)€-i ~ S 2 2(x-i,x )t; (17) 
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with £o = 1. And the unstable direction at point (xo,yo) is 
? u = (l,77 ). 

Finally, the angle between the stable and unstable directions, 
i.e. the characteristic angle of the hyperbolic structure, can 
be calculated easily 



a = arccos 



(18) 



As examples of the above algorithm, we calculated the posi- 
tions of a series of hyperbolic periodic orbits and their charac- 
teristic angles. Generally, the more irrational a rotation num- 
ber is, the more robust the corresponding KAM torus is, and 
the stronger the stickiness effect is around the relics of the 
KAM torus after its breaking. People often write the rotation 
number in the form of continued fraction, on account of the 
merit of its direct relation to the "irrationality". A continued 
fraction 

1 

aa + ; 

- B3+- 

can be denoted by [ao; a\,ci2, «3, ■ • • ], where ao is the integer 
part of the number and a, are the positive integers. Any ratio- 
nal number can be expressed as a truncated continued frac- 
tion [ao; a\ , a,i, «3, • • • , a„]. And the integer n will be called 
the "order" of the continued fraction in this paper. Note an 
« th order continued fraction ended with a„ — 2 is the same 
as an (« + 1) th order one ended with a„ - \,a n+ \ = 1. 
Several hyperbolic periodic orbits are calculated and the re- 
sults are listed in Table \T\ An obvious tendency derived 
from Table Q] is that, the bigger period number (the integer 
p of the rotation number q/p) leads to the smaller charac- 
teristic angle. As the rotation number is getting more and 
more irrational the angle becomes smaller and smaller. It 
is well-known nowadays that the last KAM tori are those 
with the "noble" rotation numbers [34], i.e. the continued 
fractions that have a* = 1 for all k above a certain num- 
ber N. And the torus with the "golden rate" rotation num- 
ber ( V5 - l)/2 = [0; 1, 1, 1, ■ ■ ■ , 1, ■ ■ ■ ] is surely the most 
robust one. In our model, we see in Table Q] that the angle 
becomes smaller when the rotation number approaches this 
noble value. And also, we notice that the "difficult-to-cross 
partial barrier" and the most sticky region we mentioned in 
§ 2.2 are just the region corresponding to this noble rotation 
number. 



4.1 Typical diffusion route 

Due to the chaotic character, a chaotic orbit in its diffusion 
route may diffuse toward any direction at a moment. Since 
the upper part of the phase space is not open to a chaotic or- 
bit but barred by KAM tori as we have seen in Fig.Q] an orbit 
will diffuse downward overall. To monitor the diffusion, we 
check the minimal y value that an orbit attains on its trajec- 
tory in the phase space. During the diffusion, each time when 
the y value reaches a new minimum y m i n , we record the time 
t. In such (f , 3> m in) records, the "back and forth" motion in the 
diffusion is abandoned, and the "one-way" diffusion is em- 
phasized. 

Arbitrarily, two initial points (xo,yo) = (0.49,29.0) and 
(0.50,30.00001) are selected and their diffusions are fol- 
lowed. The temporal evolutions of the minimums of the y 
value are shown in Fig. [7] Note only part of the orbits, when 
the minimums are smaller than y = 20.5, is displayed. Each 
point (f,y m in) m Fig. UJ tells that at time t the orbit attains a 
new minimal value y m i n . 

Rotation Number co 

2.62 2.64 2.66 2.68 2.70 2.72 




Log (Time) 

Figure 7 The variation of minimal y of an orbit with respect to time (curves 
from upper left to lower right). The red solid squares are for the orbit ini- 
tialized at (0.49, 29.0) and the black open circles for (0.50, 30.00001). The 
rotation numbers (indicated on the top x-axis) of orbits around the corre- 
sponding local lowest position are plotted (two curves crossing the panel 
from the lower left to upper right). Three regions that will be analyzed be- 
low are indicated by horizontal lines at y m j n = 20.05, 18.92 and 18.68 (see 
text for explanations). 



4 Stickiness effect and characteristic angle of 
hyperbolic structure 

When the trajectory of an orbit is plotted on the phase space, 
we can see the trajectory points accumulate in some spe- 
cific areas, indicating the existence of some structures caus- 
ing stickiness effect. In this way we may obtain a direct im- 
pression of the sticky region, just as we have shown in Fig. [3] 
Considering a bunch of numerous orbits in the phase space, 
they may leave a certain region as time passes by. The sur- 
vival probability of orbits in the definite region decays with 
time slowly. The slow decay is also an indicator of the stick- 
iness effect [16, 17, 25]. Following these techniques, in this 
section we will show the diffusion routes of some typical or- 
bits and find out where they are "stuck". We will check the 
structures causing the stickiness effect in detail, and figure 
out the relation between the diffusion speed (or inversely, the 
stickiness effect) and the characteristic angle of hyperbolic 
structures in the sticky region. 



On the (f,ymin) curve in Fig. [7] an orbit makes some steep 
"steps" and flat "plateaus". A steep step indicates that the 
orbit diffuses very quickly in this stage, probably crossing a 
well-developed chaotic region or jumping over a wide island- 
chain. While a plateau indicates that the orbit spends a long 
time achieving a small advance in y direction, i.e., it is facing 
an obstacle on its diffusion route. And the wider the plateau 
is, the more effective the obstacle is (the stronger stickiness 
effect it has). 

Owning to the chaotic character, two orbits follow different 
diffusion routes and may suffer different stickiness effects, 
that's why in Fig. [7] two curves differ from each other. But 
in the same phase space, the sticky regions that they have 
to pass through are the same, thus the plateaus appear at the 
same y m j n on two curves. Note the logarithmic scale of the 
abscissa affects the apparent widths of the plateaus. 
When an orbit reaches a new minimum y mm at moment t dur- 
ing its diffusion, a temporary (local) rotation number a> of 
the orbit around this point is calculated. For each ymin, we 
use 1001 orbital points (including the one right at moment 
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Table 1 Angles between the stable and unstable manifolds of hyperbolic periodic orbits. The first 3 columns are the rotation numbers given in different forms, 
the forth column is the order of the continued fraction, and the last column is the angle (in degrees). 



Rntntinn Mi 1 inhpr 

IMJLd.LlUll 1> LL111UC1 


pfa ptinn 
I 1 dLLUJll 


Continued Fraction 


Order 


An^le 


2.5 


5/2 


[2 


2] 


i 


169.573 


2.6666- ■ ■ 


8/3 


[2 


1,2] 


2 


116.235 


2.6 


13/5 


[2 


1.1,2] 


3 


9.034 


2.625 


21/8 


[2 


1,1,1,2] 


4 


2.094 


2.6153- • ■ 


34/13 


[2 


1,1,1,1,2] 


5 


1.042 


2.6190- • ■ 


55/21 


[2 


1,1,1,1,1,2] 


6 


1.657 


2.6176- •■ 


89/34 


[2 


1,1,1,1,1,1,2] 


7 


0.742 


2.6186- ■■ 


144/55 


r? 


1 1 1 1 1 1 1 71 


8 





f, 500 before and 500 after f) to estimate to. These re 
(oJ,ymin) are plotted in Fig. [7] too, as the curves crossir 
panel from lower left to upper right. On these curves, w 
find the temporary rotation number (w, indicated on tl 
x axis) for each y m ; n . Again, the jumping of at implie 
the orbit is temporarily in the vicinity of a big island, 
the (nearly) continuous variation of a> implies that the 
is crossing a series of finely ground island-chains, or a 
torus. The consistence between the continuous_yariatior 
and the wide plateau of y m [ n is evident in Fig. [7] 

Combiningthe explorations on the diffusion routes in 
and in Fig?|7] we know that an orbit meets considerable 
tance when crosses some regions, where the (f , jmin) c 
make plateaus and the (w, y m i n ) curves vary continuous! 
choose three such "sticky regions" to investigate in del 
the following part of this paper. Three horizontal lir 
Fig. [7] show the lower bounds of these regions. Partici 
the lowest one is close to wide plateaus on the (f , y m i n ) ci 
indicating a strong stickiness effect around it. Checkir 
rotation number, we see that it is close to the golden 
To explain clearly the "sticky region", here we defiri 
sticky region labeled Region c as an example. A lower t 
y\ = 18.68, as shown in Fig. [7] is set first. For an orbit 
ing from an initial point above this lower bound, if the y 
of this orbit gets smaller than y[, for the first time duri 
diffusion, we regard the orbit as having crossed Region 
upper bound y" is selected as well. When the y value 
orbit gets smaller than y" for the first time, the orbit is de 
as having entered Region c. The difference between th 
per and lower bound 6y c = y" - y{. is the width of this i 
region. 

It is worth to note that there is no clear/sharp boundary for 
a sticky region in the above definition. Although we use a 
Jmin value to denote it, we should keep in mind that a sticky 
region is not a straight (but most probably curved) belt across 
the phase space. 



4.2 Characteristic angle of hyperbolic structure and dif- 
fusion speed 



To explore the stickiness effects in these sticky regions, first 
we plot in Fig. [8] part of the phase space in detail around the 
corresponding local minimum y lmn where orbits feel the con- 
siderable resistance along their diffusion routes. 
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Figure 8 Three sticky regions. Three panels are for the regions indicated 
by arrows in Fig. [jj (see text). Horizontal lines in each panel represent the 
line of y = 20.05, y = 18.92 and y = 18.68, respectively. The island chains 
in red are the ones with the highest period number found in the region. 



According to Birkhoff theorem [see e.g. 31], the elliptic pe- 
riodic orbits and hyperbolic periodic orbits emerge in pairs 
after the breaking of the invariant tori in a 2D mapping. The 
elliptic periodic orbits locate in the centers of islands, while 
the hyperbolic periodic orbits exist between every two adja- 
cent islands in an island-chain. These two types of periodic 
orbits share the same rotation number. The diffusion of an 
orbit is through these hyperbolic structures. 

The elliptic periodic orbits are surrounded by invariant curves 
and they together make islands. The hyperbolic periodic or- 
bits however are unstable and the geometric structures around 
them are not as visible as the islands. So we plot in Fig.[8]as 
many as possible the island chains to imply the positions of 
hyperbolic periodic orbits and to show the composition of the 
sticky regions. In each region, the island chain with the high- 
est periods (they are 71, 73 and 144 in Region a, b and c) 
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Table 2 Characteristic angles of hyperbolic structures in three sticky regions. The hyperbolic periodic orbits are denoted by their rotation numbers, given in 
normal fraction numbers and continued fractions. The angles are in degrees. 



Region a 



Rotation number 


Angle (a) 


Fraction 


Continued Fraction 


165/61 


[2;1,2,2,1,1,3] 


0.718653 


119/44 


[2; 1,2,2, 1,1, 2] 


0.845155 


192/71 


[2; 1,2,2, 1,1, 1,2] 


0.772372 


73/27 


[2;1,2,2,1,2] 


1.094884 


Region b 


171/65 


[2;1, 1,1,2,2,3] 


0.717274 


121/46 


[2; 1,1, 1,2,2,2] 


0.863938 


192/73 


[2;1, 1,1,2,2,1,2] 


0.777229 


71/27 


[2;1, 1,1,2,3] 


1.179671 


Region c 


199/76 


[2;1, 1,1, 1,1, 1,1,3] 


0.506420 


343/131 


[2;1, 1,1, 1,1, 1,1,2,2] 


0.472460 


144/55 


[2;1, 1,1, 1,1, 1,1,2] 


0.582802 


377/144 


[2;1, 1,1, 1,1, 1,1, 1,1,2] 


0.477654 


233/89 


[2;1, 1,1, 1,1, 1,1, 1,2] 


0.534822 


89/34 


[2;1, 1,1, 1,1, 1,2] 


0.741975 



have been plotted in red. Of course, there must be islands of 
even higher periods, but their sizes are much smaller, and it 
is impossible to plot all of them in practice. The hyperbolic 
periodic points can be located and the characteristic angles 
of the hyperbolic structures can be calculated, using the algo- 
rithm introduced in above section. Some characteristic angles 
in these regions are summarized in Table [2] 
We have seen from Table Q] that the characteristic angles of 
hyperbolic structures with higher period numbers are smaller. 
This trend can be seen again in each region in Table [2] But 
when considering the characteristic angles with nearly the 
same period number in different regions, for example, the 
characteristic angles corresponding to the rotation number 
192/71 in Region a, 192/73 in b and 199/76 in c, we find that 
the angle in Region c (0.506°) is significantly smaller than 
the other two in Regions a (0.772°) and b (0.777°), while the 
latter two are nearly equal to each other. This result is ex- 
pectable since Region c makes the widest plateau in Fig. [7] 
and the diffusion time through this region is in the order of 
10 8 ~ 10 9 , much longer than that in Regions a and b. 
The characteristic angle may decreases as the period of orbits 
increases, but it is impossible in practice to compute the char- 
acteristic angles for all the hyperbolic structures. Fortunately, 
the information about the smallest characteristic angle in a re- 
gion can still be obtained. In Region c, the rotation numbers 
of the hyperbolic periodic orbits is close to the "golden rate", 
and it is well-known that the KAM torus with the golden rate 
rotation number is the most robust one. Meanwhile, it is nat- 
ural to assume that the "smallest" characteristic angle should 
be found in the hyperbolic structures of the new-born can- 
torus just after the breaking of KAM torus. Since the golden 
rate can be approached by continued fractions, the character- 
istic angles of hyperbolic structures with rotation numbers of 
truncated continued fractions may approach to the character- 
istic angle in the cantorus. Following this idea, we continue 
the calculations listed in Table Q] to higher order of the con- 
tinued fraction. The variation of the characteristic angle with 



respect to the order of the continued fraction is illustrated in 
Fig. [9] Note, as the order increases the corresponding peri- 
odic orbits tend to gather around the golden rate cantori, and 
no later than the 7th order (89/34) all the hyperbolic struc- 
tures under consideration accumulate in Region c. 




Order of Continuous Fraction 

Figure 9 The lower bound of characteristic angles in Region c. The ab- 
scissa is the order of the truncated continued fraction, which is the rotation 
number of the hyperbolic periodic orbit. The ordinate is the corresponding 
characteristic angle. The last point of order 16 represents the hyperbolic 
periodic orbit with rotation number 10946/4181. 

As shown in Fig. [9] the calculations are stopped at the 16th 
order where the rotation number is 10946/4181. We cease 
the calculations for two reasons as follows. First, the calcu- 
lations so far have shown that the characteristic angle has a 
lower bound when the truncated continued fractions approach 
the golden rate. Even higher order does not decrease this 
lower bound. Second, calculations for the high order contin- 
ued fractions (therefore high period number) are very expen - 
sive in computer time, and it is difficult to control the error 
in this case because we have to deal with a large number of 
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linear equations (Eq. (15)). So we adopt the smallest angle in 
Fig- H] at the 12th order, a c = 0.457°, as the lower bound of 
the characteristic angle in Region c. 

For Regions a and b, there is no golden rate rotation num- 
ber. However we can follow a similar process as we have 
done for Region c to calculate the lower bound of character- 
istic angles. We expand the continued fractions listed in Ta- 
ble |2]by inserting number '1' into them right before the last 
number, then we have series of truncated continued fractions 
of different orders. Adopting these continued fractions as the 
rotation numbers, we calculate the characteristic angles of the 
corresponding hyperbolic structures. In this way, we obtain 
the lower bound of the characteristic angles a a = 0.665° in 
Region a and a/, = 0.667° in Region b. 
Since a a « a/, > a c , we may derive directly from these values 
that the diffusion in Region c is slower than that in Regions 
a and b. In other words, the stickiness effect in Region c is 
more intensive than that in Regions a and b, while Regions a 
and b have nearly the same stickiness effects. To verify this 
conclusion, we carry out some calculations of diffusion speed 
below. 




Time Tim e 

Figure 10 The diffusion speeds of orbits in Regions a, b and c. The time 
is given in linear scale in left panel, and logarithmic scale in right panel. 

The diffusion speed, generally is measured by the decay rate 
of the survival numbers of a bunch of orbits in a defined area 
[8, 16, 17]. We set randomly 40,000 initial points in a rect- 
angle from (0.45, 26.5) to (0.55, 31.5) around the hyperbolic 
fixed point of rotation number 3/1 that is embedded in a de- 
veloped chaotic area and above the Regions a, b and c in 
the phase space (see Fig. []]). All the orbits are iterated for 
10 9 times using the mapping. The evolution of these forty 
thousand orbits are then followed, especially, we monitor the 
variation of the action variable y of each orbit. 
The definition of sticky region has been introduced in § 4. 1 . 
For Regions a, b and c, the lower boundaries are chosen to be 
y' a = 20.05, y' b = 18.92, y' c = 18.68 respectively, as indicated 
by the horizontal lines in Figs.|7]and|8] After some tests, we 
set the width of the sticky region 5y = 0.1 for all three re- 
gions. The upper boundaries of the sticky regions are then 
defined as y" = y l , + 5y. 

An orbit may wander for quite a while before it enters a spe- 
cific area. To exclude the evolution history before an orbit en- 
tering the vicinity of Region a, b or c, we define the moment 
f a (or t" b , t") when for the first time the y value of an orbit is 
smaller than the upper boundary y" (or y b ,y") as the moment 
of entering Region a (or Region b, c). On the other hand, 
when an orbit satisfies for the first time y < y l a (or y' b ,y l c ) at 
moment t' a (or H , f'), it is regarded as having escaped from the 
sticky region. The time duration At a ^,c = K,hc~^abc tnen * s 
defined as the surviving time of an orbit in the corresponding 
sticky region. Following the evolution of the 40,000 orbits, 



we record the surviving time of each orbit in three sticky re- 
gions. With these records, we count how many orbits survive 
in a specific region at a given surviving time. Some orbits 
may have been initialized on invariant curves by chance, and 
a small number of orbits may wander in the chaotic phase 
space but never reach a specific region. These orbits are ex- 
cluded from the final statistics. After removing these orbits 
and normalizing the number of orbits at the beginning to unit, 
the diffusions in the three regions are summarized in Fig. HQ] 
Apparently, the diffusion is quick in Regions a and b, but 
slow in Region c, reflecting that the stickiness effect in Re- 
gion c is much stronger than in other regions. On the other 
hand, although Regions a and b occupy different areas in 
the phase space and they comprise different components as 
shown in Fig. [8] they have nearly the same diffusion speed. 
The results are consistent with the conclusion we derived 
from the characteristic angle values. 

4.3 Characteristic angle and perturbation parameter 

So far we have seen that in the phase space of a mapping with 
given perturbation parameter k, the stickiness effects of dif- 
ferent sticky regions are related to the characteristic angles of 
hyperbolic structures embedded in the regions. Smaller char- 
acteristic angles lead to stronger stickiness effects and slower 
diffusions. Below we will investigate the relation between 
the stickiness effect and the characteristic angle when the pa- 
rameter k changes. We will focus on the Region c where the 
diffusing orbits suffer the most significant stickiness effects 
as shown above. 




/ 19 20 21 22 23 

Parameter (£) 

I i I i I i I i I i I i 1_ 

20 25 30 35 40 45 50 

Parameter (k) 



Figure 11 The variation of characteristic angle with respect to the param- 
eter k. The solid circles show the angle of the hyperbolic structure with ro- 
tation number 199/76. The embedded panel shows the characteristic angles 
of several hyperbolic structures with different rotation numbers. 

In Region c, we select five hyperbolic periodic orbits to an- 
alyze, including the one with rotation number 377/144 for 
which the corresponding island-chain was plotted in red in 
Fig. [8] three above it with rotation numbers 199/76, 343/131 
and 144/55, and one below it with the rotation number 
233/89. The corresponding island chains are all plotted in 
Fig. [8] and they stay close to each other in a narrow area in 
the phase space. Their characteristic angles at different k val- 
ues are then calculated . The results are illustrated in Fig. QT| 
and some of them are listed in Table [3] too. 
When k is small, the global KAM ton exist in Region c and 
these hyperbolic structures are bounded by these tori. As 
pointed out in several literatures [e.g. 27, 36], in the vicin- 
ity of an invariant torus in a 2D phase space, the resonances 
(island-chains and hyperbolic periodic orbits) accumulate ge- 
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Table 3 The characteristic angles at different parameter k. The angles are in degrees. Only five hyperbolic structures in Region c are listed here. 



Rot. Num. 


Angle (a) 


k = 19.00 


20.35 


21.50 


22.50 


27.50 


32.50 


36.00 


199/76 


0.396155 


0.448865 


0.792791 


1.291957 


3.288836 


11.42482 


77.06509 


343/131 


0.429506 


0.426083 


0.799331 


1.294550 


3.288836 


11.42482 


77.06509 


144/55 


0.405655 


0.504247 


0.851525 


1.295668 


3.288836 


11.42482 


77.06509 


377/144 


0.431027 


0.437425 


0.840257 


1.314751 


3.288844 


11.42482 


77.06509 


233/89 


0.431263 


0.467751 


0.856823 


1.314805 


3.288844 


11.42482 


77.06509 



ometrically toward the torus with their increasing orders (pe- 
riod numbers). Owning to continuity, the directions of the 
stable and unstable manifolds of the hyperbolic structures are 
"forced" to be nearly parallel to the torus, thus the charac- 
teristic angles are quite small in the close neighborhood of a 
torus. Here we see from the embedded panel in Fig. [TTIthat 
the angles seem to approach a limit slowly as k decreases. 
Consequently, orbits starting from the close vicinity of an in- 
variant torus will diffuse along the resonance line (high-order 
island-chain) for a long time before its leaving for adjoining 
resonances. As k increases, the chaos is developed and the 
characteristic angle increases. Particularly, after the global 
KAM tori disappear in this region (when k > 20.75) the char- 
acteristic angle increases dramatically. As a result, the diffu- 
sion speed is expected to be enhanced dramatically. 
Another interesting phenomenon is the difference between 
the characteristic angles of different rotation numbers. When 
k is small, the characteristic angles are small and the rela- 
tive difference between different rotation numbers is consid- 
erable. But the relative difference becomes smaller as the 
angles increase and almost invisible when a ~ 2°(k = 25). 
This is clearly illustrated in the embedded picture in Fig. fTTl 
The reason is that the difference of hyperbolicity between dif- 
ferent periodic orbits is apparent when k is small, but it be- 
comes ignorable when k is large. In the latter case, chaos has 
developed, making the whole region uniformly hyperbolic. 
Thanks to this observation, to show the angle variation in a 
wide k range, it is enough to follow just one of the hyper- 
bolic structures, e.g., the one shown in Fig. QT|with rotation 
number 199/76. 



4.4 Stickiness effect and characteristic angle 

The stickiness effect (or the diffusion speed) changes when 
the characteristic angle changes with k. We will present the 
relation between the diffusion speed and the size of charac- 
teristic angle in this part. To show this correspondence, we 
need to define the diffusion speed at first. The time needed 
by a bunch of orbits to cross through a definite sticky region 
may be a good measure. But a "sticky region" may deform as 
the structure of phase space distorts with the varying param- 
eter k. Thus we check first the variation of the position of a 
hyperbolic periodic orbit at different k. 
We plot in Fig. Q~2] the hyperbolic periodic orbit of rotation 
number 377/144 at six k values. As k increases, the profile of 
the hyperbolic periodic orbit is distorted. The other four hy- 
perbolic periodic orbits of rotation numbers 199/76, 343/131, 
144/55 and 233/89 are calculated too, and we find that their 
profiles are distorted in a similar way as shown in Fig. Upland 
they always stay together within a narrow area in the phase 
space. Therefore, the "sticky region" for different k can be 
defined in a similar way as in § 4.2. 



if 
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Figure 12 The positions of the hyperbolic periodic orbits of rotation num- 
ber 377/144. Cases of different k are represented by different symbols. 

Just as we have done for the case k = 20.75, forty thousand 
points initialized in the rectangle area are followed again for 
different k values. We now focus only on the Region c where 
the above mentioned five hyperbolic periodic orbits are em- 
bedded. Knowing the locations of these hyperbolic structures 
(as shown in Fig. [12b, the "sticky region" at different k is de- 
fined again. And also, the surviving time of an orbit in the 
sticky region is given by the time duration At - t l - f as 
before. 




Figure 13 The diffusion at different k. From top to down (right to left), the 
curves represent the case for k = 20.75, 21.00, 21.50, 22.50, 25.00, 27.00 
and k = 30.00, respectively. 

All the orbits are iterated for 10 9 times. We record the dura- 
tion time At of each orbit, and it is regarded as the surviving 
time of the orbit in the sticky region. The number of surviving 
orbits changes with time. Following the same technique for 
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plotting Fig. [TOb . we summarize our calculations for several 
k in Fig. |T3l 

It is clear that the diffusion speeds up as k increases. For 
the case of k = 20.75, the first orbit escaping happens at 
t ~ 2 x 10 5 , when more than 90 percent of orbits have es- 
caped in the case of k — 22.50. Nearly all orbits of k — 25.0 
and k = 27.00 escape before / ~ 10 4 , while only 1 percent of 
orbits survive after only hundreds of iterations for k = 30.00. 
In many literatures [e.g. 8, 16, 17, 35], the diffusion curves 
similar to the ones in Fig. [13] were fitted by a power law 
/^survival oc t~ z . The exponent z then can serve as a measure- 
ment of the diffusion speed (or the intensity of the stickiness 
effect). But such a power law is not so distinct in our cases as 
we see in Fig. [13] Instead, we use directly the time when most 
of the orbits have escaped to measure the diffusion speed, (or 
inversely, the intensity of stickiness effect). More or less arbi- 
trarily, we adopt the moment when 90 percent of orbits have 
escaped, denoted by t\/io since only one tenth of orbits sur- 
vive at this moment, as an estimation of the diffusion speed. 
Thus a quick diffusion has a small fi/io while a large t\/\o 
indicates a slow diffusion (and strong stickiness effect). 
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Figure 14 The variation of escaping time with respect to the characteristic 
angle (a). The abscissa is cot a = 1/tana and the ordinate is the escap- 
ing time (fi/io, see text for definition) in logarithm. Several a values are 
indicated along the top axis too. 

For k varying from 20.75 to 30.00, we calculate the ti/io and 
the characteristic angles in this region at these k values are 
computed as well. They are illustrated in Fig. [14] A small 
characteristic angle leads to a slow diffusion, that is, a large 
fi/io- So, to visualize the variation of diffusion speed at small 
characteristic angle, we use 1 / tan a instead of a itself as the 
abscissa in Fig. [14] while several a values are indicated along 
the top axis. The ordinate is the diffusion speed t\/io in loga- 
rithm scale. For a given parameter k, the characteristic angle 
in Fig. [14] is the lower bound of the characteristic angles in 
Region c. They are calculated in the same way that we have 
introduced in § 4.2. 

As indicated by Fig. [14] the stickiness effect is significant 
only when the characteristic angle is small. When the char- 
acteristic angle a ~ 5°(1/ tana- = 11.5), most of orbits have 
escaped just after tens of iterations, implying a weak stick- 
iness effect. On the other end of smaller a in Fig. [T4J the 
escaping time will be longer than 10 9 when a is smaller than 
0.45°(1/ tan 0.45° = 127.3). An orbit with such a long dif- 
fusion time can be practically regarded as stable (never es- 
cape). In fact, the global KAM curves exist in Region c when 
k < 20.5. And the characteristic angle is smaller than 0.45° 
when k < 20.5, as we can see from Fig.fTTI 



It is worth to note that the points in Fig. [14] may be fitted by 
two linear functions with a crossover around a ~ 1°. The 
possible mechanism underneath such a relation deserves an 
investigation in future. 

When k > 20.75, all the global KAM tori in Region c have 
broken. The stickiness effect suffered by an orbit crossing 
this region arises from the island-chains and cantori embed- 
ded in the region. The diffusion speed is then determined by 
the characteristic angle of hyperbolic structures among the 
island-chains and cantori. Generally, in the vicinity of invari- 
ant tori, hyperbolic structures accumulate, thus the stickiness 
effect felt by orbits close to the tori, where the stickiness ef- 
fect was firstly found and defined [?, ?], has the same origin: 
the hyperbolic structures. 

5 Conclusions 

The stickiness effect in the phase space is an interesting phe- 
nomenon with implications in many areas of sciences. It has 
attracted many attention in the fields of mathematics, physics 
and astronomy. Many structures such as the KAM tori, the 
island chains, the cantori, and the hyperbolic structures in the 
phase space are found to have stickiness effect. Among them, 
we think the hyperbolic structures play the essential role, as 
we have shown in our previous papers. 
In this paper, we describe some geometric details of the hy- 
perbolic structure, especially the angle between the stable 
and unstable manifolds, and relate these geometric characters 
with the strength of the stickiness effect. 
Usually the recursive algorithm is employed to calculate the 
positions of periodic orbits in the phase space. But a usual 
recursive algorithm fails when we try to find the hyperbolic 
periodic orbits, because the set of hyperbolic periodic orbits 
is of measure zero when the system is far from the integrable 
one. In this paper, a numerical algorithm is introduced, with 
which we can compute the precise locations of the hyperbolic 
periodic orbits. And also, a practical algorithm for computing 
the directions of stable and unstable manifolds of the hyper- 
bolic structures is presented. 

With these numerical algorithms, we have computed the an- 
gles between the stable and unstable manifolds in different 
hyperbolic structures in the 2D phase space of a mapping 
model. We investigate the stickiness effect (diffusion speed) 
in different regions of the phase space, and how it changes 
with the perturbation parameter. 

Our findings in these calculations may be summarized as be- 
low: 

- In an area where the last KAM torus has broken, the char- 
acteristic angle between the stable and unstable manifolds of 
the hyperbolic structure is related to the rotation number. The 
more irrational the rotation number is, the smaller the angle 
is. But a lower bound of this angle exists in a certain region. 

- The stickiness effect (diffusion speed) is determined by the 
characteristic angle. The smaller the angle is, the slower the 
diffusion is, alternatively speaking, the stronger the stickiness 
effect it has. 

- A relationship between the characteristic angle of the hy- 
perbolic structure and the stickiness effect has been revealed. 
Nevertheless, it is still difficult to find a quantitative way to 
describe this relationship explicitly. 

Although the above conclusions are drawn from the calcula- 
tions in regions consisted of island-chains and cantori but not 
global invariant tori, they should be valid for all structures 
proven to possess stickiness effect. 
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